Superconducting Proximity Effect and 
Universal Conductance Fluctuations 
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We examine universal conductance fluctuations (UCFs) 
in mesoscopic normal-superconducting-normal (N-S-N) struc- 
tures using a numerical solution of the Bogoliubov-de Gennes 
equation. We discuss two cases depending on the presence 
("open" structure) or absence ("closed" structure) of quasi- 
particle transmission. In contrast to N-S structures, where 
the onset of superconductivity increases fluctuations, we flnd 
that UCFs are suppressed by superconductivity for N-S-N 
structures. We demonstrate that the fluctuations in "open" 
and "closed" structures exhibit distinct responses to an ap- 
plied magnetic field and to an imposed phase variation of the 
superconducting order parameter. 



Transport properties of mesoscopic structures in con- 
tact with superconducting segments have attracted a 
great deal of interest In these structures, the 

superconducting proximity effect induces novel trans- 
port phenomena, including sub-gap anomalies in normal- 
superconducting (N-S) contacts and phase-dependent 
conductances in Andreev interferometers. Quasiclassical 
theory, which has been widely applied, has been success- 
ful in explaining phenomena associated with ensemble- 
averaged properties but it cannot describe fluctuations 
about the mean 

In normal diffusive nanostructures of size smaller than 
the phase coherence length but larger than the elas- 
tic mean free path l^i, the conductance fluctuates by an 
amount of order e'^/h. These fluctuations can be ob- 
served by varying either the magnetic field, gate poten- 
tial, or impurity configuration. The fluctuations are uni- 
versal within the diffusive regime: their magnitude is 
partly determined by symmetry and does not depend on 
the degree of disorder or system size. Such a symmetry- 
breaking operation is the application of an external mag- 
netic field which breaks time reversal (TR) symmetry and 
decreases the rms fluctuations by a factor of V2. 

For normal-superconducting systems, it has been 
shown that the onset of superconductivity increases 
the rms magnitude of UCFs by a factor which depends on 
the underlying order parameter. Beenakker and Brouwer 
(BB) 1^ predict that in the absence of phase gradients for 
the order parameter, breaking TR symmetry has a negli- 



gible effect on the rms conductance fluctuations (5Gns of 
a N-S structure. Altland and Zirnbauer (AZ) on the 
other hand, find that (JGns decreases by a factor of \/2, 
provided that the average phase change of an Andreev- 
reflected quasiparticle is zero. When the latter condition 
is satisfied, AZ predict that in zero magnetic field. 



SGm{B = 0) = 2\/2<5Gnn(S = 0), 



(1) 



where ^Gnn is the rms conductance fluctuation when 
the S contact is in the normal state. Since switching on 
a magnetic field decreases both 5Gns and 5Gnn by a 
factor of v^, they also find 

(5Gns(S 7^ 0) = 2V2JGnn(S ^ 0) = 2(5Gnn(S = 0). 

(2) 

In contrast, since BB predict that (5Gns is almost insen- 
sitive to TR symmetry breaking, they obtain 

SG^siB = 0) « 26GMB = 0) (3) 

SG^siB 7^ 0) = 2%/2^Gnn(S 7^ 0) = 2(5Gnn(S 0). 

(4) 

Thus the AZ and BB scenarios yield essentially the same 
fluctuations for finite B, a result which has recently been 
tested experimentally |g] , whereas they yield fluctuations 
of different magnitudes in zero field. 

In this paper we examine for the first time the crossover 
between the AZ and BB scenarios and predict that this 
crossover is observable in closed Andreev interferometers. 
The BB and AZ calculations are applicable to N-S struc- 
tures with a single N contact, whereas many experiments 
involve two or more N contacts Q . We also examine the 
field and phase dependence of fluctuations in the conduc- 
tance Gnsn of N-S-N structures possessing two normal 
(N) contacts [0,^. For the considered geometries, these 
have not been explored before in the literature. 

We consider the two limiting cases of N-S-N structures 
shown in Fig. 1, comprising two N-reservoirs attached to 
a normal diffusive region with superconducting contacts 
and a current flowing from left to right. In Fig. la, the 
two superconductors do not intersect the diffusive nor- 
mal region and affect transport only through the prox- 
imity effect. In this structure, quasiparticle transmission 
between the reservoirs is possible and therefore we refer 
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to this as "open". In the second, "closed", structure of 
Fig. lb, a superconductor of length Ls divides the diffu- 
sive region into two normal regions, each of length L^. 
The length Lg of the superconducting segment is much 
greater than the superconducting coherence length ^ and 
therefore quasiparticle transmission through the super- 
conductor is suppressed. 






FIG. 1. a) An open Andreev interferometer where the su- 
perconducting segments are located outside the classical path 
of the current. The N-S interfaces contain an optional Schot- 
tky barrier with transmissivity F (i.e., transmission probabil- 
ity per channel). The S segments have an order-parameter 
phase difference of = <j!>i — <j!>2 b) A closed interferometer 
where the superconducting segments divide the normal dif- 
fusive region into two halves with no transmission from left 
to right. At the interface between the S segments there is a 
tunnel barrier of width 1 site and transmissivity F = 0. 

In the linear-response limit and at zero temperature, 
the relationship between the current / and the reservoir 
potential difference vi — V2 was first obtained in Q and 
the corresponding conductance (in units of 2e^/h) may 
be written in the form [|1|, 



G 



NSN 



= To + Ta + 2 



Ra+R'+Ta + T' 



(5) 



where unprimed [primed] quantities Tq and i?o (T^ and 
Ra) are normal (Andreev) transmission and reflection 
probabilities for quasiparticles from lead 1 [lead 2]. To 
evaluate this expression, we solve the Bogoliubov-de 
Gennes equation numerically on a two-dimensional tight- 
binding lattice with diagonal disorder, utilizing the dec- 
imation method to evaluate the transport coefficients of 
a two-probe structure The magnetic field is intro- 

duced via the Peierls substitution and the disorder in the 
normal diffusive area is introduced by varying the site 
energies at random within the range [— it;/2, ■u;/2]. 

Conductance fluctuations are universal only if trans- 
port through the normal structures is diffusive. There- 
fore, following [0, we choose the disorder w such that 
the conductivity of the diffusive normal region is inde- 
pendent of its length L. For a sample of width M = 35, 
in the absence of superconductivity. Fig. 2 shows the di- 
mensionless conductivity I = {G)L/M as a function of 
L for four different values of disorder w. In the quasi- 
ballistic regime, (G) approaches the Sharvin conductance 
and I increases linearly with L, whereas at intermediate 
L, / exhibits a plateau, corresponding to the diffusive 



regime of interest. In units of the lattice constant, this 
plateau value of I is equal to the elastic mean free path 
Ici- Therefore, Fig. 2 represents a mapping between the 
model-dependent parameter w and the experimentally 
accessible quantity Id . Finally, for values of L larger than 
the localization length A « Nlc\, where = 30 is the 
number of open channels, I decreases with increasing L. 
In Figs. 4-5 we choose w = 1.5, corresponding to a mean 
free path of ^oi ~ 8. The other length scales (in units of 
the lattice constant) are M = 35, Ln = 75, Ls = 40, 
W — 10 and — 10, where ^ is the superconducting 
coherence length setting the scale for the decay of zero- 
energy wavefunctions within the superconductor . 

20- 
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FIG. 2. Quantity {G)L/Ai plotted as function of L for 
different strengths w of disorder. The region where the curve 
is essentially constant corresponds to diffusive transport. The 
constant defines the mean free path l^i- The dotted lines at 
I = Id are guides to the eye. 

Figure 3 shows the rms fluctuations as a function of 
I for a normal structure in the absence of superconduc- 
tivity and for the N-S-N structures of Fig. 1 (with zero 
magnetic field and a constant order-parameter phase) ; it 
demonstrates that fluctuations are independent of This 
universal behavior is in agreement with previous numer- 
ical and analytical [^ Q works. In the diffusive 
regime, the presence of superconductivity for both struc- 
tures in Fig. 1 suppresses the fluctuations by a factor of 
order \/2- 

For the closed structure in Fig. lb, all transmis- 
sion coefficients necessarily vanish, and for structures 
of type la, the presence of disorder and a proximity- 
induced suppression of the density of states may also 
cause transmission coefficients to be negligible. In this 
case, Eq. (H) reduces to a sum of two N-S resistances, 
Gnsn = yC^Ra) + l/(2-Ra), associated with Andreev 
reflection of quasiparticles from the two decoupled reser- 
voirs. By adding two statistically independent N-S resis- 
tances in series, one flnds that in the absence of transmis- 
sion, the fluctuations (5Gnsn of the total conductance are 
related to the fluctuations SG^s of the N-S-conductance 
by (5Gnsn = i5Gns/2'\/2. Therefore, in the absence of 
transmission, the decrease in i5Gnsn due to the onset 
of superconductivity is a consequence of this relation and 
Eq. (3) II . 
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FIG. 3. Universal conductance fluctuations (rms devia- 
tions) in units of 2e^//i as function of I for a normal structure 
(crosses), open N-S-N structure (open circles), and closed N- 
S-N structure (closed circles), obtained from over 500 real- 
izations of disorder. The mean free path extends over three 
regimes: localized (/ < 6), diffusive {6 < I < 35) and qua- 
siballistic (/ > 35). The solid lines at 5G = SGn = 0.358 
and SG — SGs ~ 0.267 show the average of the standard 
deviations within the diffusive regime in the normal and su- 
perconducting cases, respectively. The dashed lines indicate 
95% confidence intervals. 

Figure 4a shows that in the presence of a large enough 
magnetic field, the rms fluctuations of both a normal 
structure and the open structure in Fig. la are suppressed 
by a factor of order ^/2. However, in agreement with the 
prediction of BB, the suppression of ^Gnsn for the closed 
structure of Fig. lb is at most of the order of 10%. 

To demonstrate that a crossover from the BB to the AZ 
scenario is possible in such interferometers, we have also 
investigated the effect of imposing a difference (j) between 
the order parameter phases of the two superconductors. 
The results of these calculations are shown in Fig. 4b. 
For the closed structure, in the absence of magnetic field 
(solid line) , (JGnsn increases monotonically as increases 
from zero, reaching a maximum at = tt. The enhance- 
ment factor at (/) = TT is of order v^, which is the ratio 
of the AZ and BB predictions for UCFs, thereby demon- 
strating that a crossover is indeed possible in closed inter- 
ferometers. In the presence of a magnetic field (dashed 
line), the fluctuations become nearly independent of the 
phase difference (fi since time reversal symmetry is already 
broken by the field. Remarkably, the corresponding curve 
for the open structure shows no phase dependence. At 
present, there exists no analytic derivation of this result. 
In contrast with the mean conductance (G), which at 
i? = is independent of 0, but exhibits large-scale os- 
cillations at finite E [||,|l6|, we have found no significant 
phase dependence in the fluctuations, even at energies of 
the order of the Thouless energy. 

Figures 3 and 4 demonstrate that if the N-S interfaces 
are clean enough, the magnitude of UCFs are sensitive 
to the onset of superconductivity in the S contacts. As a 
final result, we note that this effect is suppressed by the 
presence of Schottky barriers at the N-S interfaces. This 
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FIG. 4. a) External magnetic field breaks time reversal 
symmetry and decreases the magnitude of fiuctuations by a 
factor of 1.3 both in a normal structure (crosses) and in an 
open interferometer (open circles). However, for the closed 
structure of Fig. lb (solid circles), the suppression is only of 
the order of 10%. The dashed lines at 6Gnsn{B / 0) = 0.207, 
5Gnsn(B ^ 0) = 0.254 and (5Gnn(B 0) = 0.275 show the 
average of fluctuations for fluxes $ > $o/2 = h/2e in open, 
closed, and normal structures, respectively, b) Magnitude 
i5Gnsn of universal conductance fiuctuations as functions of 
the phase difference <j> — 4>i — 4i2 between the superconduct- 
ing segments in the absence (solid line and circles) and the 
presence (dashed line and triangles) of a magnetic field. Re- 
sults are presented for both closed (solid markers) and open 
interferometers (open markers). 

is demonstrated in Fig. 5 which shows (5Gnsn as a func- 
tion of the barrier conductance Gb — NT for different 
numbers of channels N in the contact and transmissiv- 
ities r. For Gb ^ 10e^//i, the fluctuations approach 
those of a purely normal structure and for Gb ^ 20e^//i 
fiuctuations characteristic of a N-S-N structure with no 
barriers are obtained. From an experimental viewpoint, 
this suppression could be used as a novel probe of the 
strength of such barriers. 

In summary, we have used the numerical solution of 
the Bogoliubov-de Gennes equations to demonstrate that 
conductance fluctuations in normal-superconducting 
mesoscopic structures depend on (i) the presence or ab- 
sence of quasiparticle transmission through the structure, 
(ii) the external magnetic fleld, and (iii) the order-para- 
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FIG. 5. Tunnel barrier at the N-S interface of the open 
structure in Fig. la weakens the effect of superconductivity 
and for barrier conductance Gb ^ lOe^/h, the fluctuations 
start to resemble those in the absence of superconductivity. 
Here the barrier conductance is expressed in the units of 2e^//i 
and fluctuations are plotted for 47 (solid square markers), 38 
(open squares), 35 (solid circles), 23 (open triangles) and 19 
(solid triangles) open channels at the interface. 

meter phase difference between the superconducting con- 
tacts. If the phase of the superconducting order parame- 
ter is constant and time-reversal symmetry is not broken 
by a magnetic field, the open and closed structures pos- 
sess the same UCFs, namely (5Gnsn ~ 0.26[2e^//i]. In 
contrast to the predictions of AZ and BB for N-S struc- 
tures, this is smaller by roughly a factor of \/2 than the 
normal-state fluctuations (5Gnn- Applying a magnetic 
field to the open structure decreases (5Gnsn by an addi- 
tional factor of \/2, a result which lies outside current 
analytic calculations ||l^, whereas for the closed struc- 
ture the suppression is found to be much weaker (roughly 
10 %), as suggested by BB. We have also examined the 
phase dependence of UCFs and shown that for a closed 
structure, (5Gnsn increases by almost \/2 as the phase 
difference between two superconducting segments is var- 
ied from zero to tt, thereby demonstrating a crossover 



System 


5GNSN/[2eV/i] 


Normal, T 


0.358 


Normal, no T 


0.275 


Open, (/> = 0, T 


0.266 


Open, <j!> = 0, no T 


0.207 


Open, (j> — n, T 


0.264 


Open, (j) = n, no T 


0.201 


Closed, = 0, T 


0.273 


Closed, = 0, no T 


0.254 


Closed, (j) = n, T 


0.336 


Closed, (j) — -K, no T 


0.255 



TABLE I. Calculated universal conductance fluctuations 
(std(G)) for normal structures, open Andreev interferome- 
ters and two N-S contacts in series ("closed"). In the latter, 
the conductance fluctuations through a single contact are ob- 
tained by a multiplication by \/8. T denotes time-reversal 
symmetry. The relative error of the magnitudes of fluctua- 
tions is obtained from 95 % confidence intervals for standard 
deviation and is of the order of 5 %. 



from the BB scenario to the AZ scenario. In contrast 
for the open structure, we find that (5Gnsn is essentially 
independent of phase. The calculated magnitudes of the 
UCFs are summarized in Table ||. 

The numerical simulations have been performed in a 
Cray C94 of the Center for Scientific Computing (CSC, 
Finland) . TH acknowledges the postgraduate scholarship 
awarded by Helsinki University of Technology and the 
hospitality of the Lancaster University. 
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